Geometrical edge barriers and magnetization in superconducting strips with slits 
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We theoretically investigate the magnetic-field and current distributions for coplanar supercon- 
ducting strips with slits in an applied magnetic field H a . We consider ideal strips with no bulk 
pinning and calculate the hysteretic behavior of the magnetic moment m y as a function of H a due 
solely to geometrical edge barriers. We find that the m y -H a curves are strongly affected by the slits. 
In an ascending field, the m y -H a curves exhibit kink or peak structures, because the slits prevent 
penetration of magnetic flux. In a descending field, m y becomes positive, because magnetic flux is 
trapped in the slits, in contrast to the behavior of a single strip without slits, for which m y « 0. 

PACS numbers: 74.25.Op, 74.25.Qt, 74.25.Ha, 74.78.-w 
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II. COMPLEX FIELD AND MAGNETIC 
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Superconducting flat strips subjected to a perpendic- 
ular magnetic field show magnetic hysteresis, even when 
the strips have no bulk pinning. The magnetic hystere- 
sis of strips without bulk pinning arises from geometrical 
edge barriers, i.e., barriers for magnetic-flux penetration 
at the strip edges Current-carrying strips have 

finite critical currents arising from the edge barriers^ and 
the critical current becomes larger when slits are fabri- 
cated near the edges of the strips^ The critical-current 
increase is due to the enhancement of edge-barrier ef- 
fects; in other words, making slits increases the number of 
edges that prevent flux penetration into the inner strips. 
The reversible magnetic response of two strips (i.e., a 
strip with a slit) in the Meissner state was considered in 
Refs.iandE When an applied magnetic field exceeds 
a certain value, magnetic flux penetrates into the strips 
and the magnetic response becomes irreversible and hys- 
teretic m- 



In this paper, we present a theoretical investigation of 
the magnetic hysteresis of bulk-pinning-free strips with 
slits in the presence of an applied magnetic field H a . 
The hysteretic behavior of the magnetic moment m y as 
a function of H a shown in this paper is due solely to the 
geometrical edge barriers. In Sec. II we outline our the- 
oretical approach and establish notation. In Sec. IIIII we 
briefly review published results of the m y -H a curves of a 
single strip without slits. In Sec. II VI we investigate field 
distributions and m y -H a curves of two strips (i.e., a strip 
with a slit), and in Sec. |V|we study three strips (i.e., a 
strip with two slits). We briefly summarize our results in 
Sec. VI. 



We investigate coplanar superconducting strips (i.e., 
strips in which slits are fabricated parallel to the edges) 
in a perpendicular magnetic field but carrying no net 
transport current. The strips under consideration have 
total width 2a, thickness d <C 2a, and infinite length 
along the z axis (i.e., |x| < a and \y\ < d/2), as shown 
in Fig. n We assume that magnetic flux penetrates and 
escapes along the x axis, assuming that there is no flux 
penetration from the ends at \z\ — * oo. (We may think 
of the strip ends as being connected by superconducting 
shunts.) 

The Biot-Savart law for the complex fielc&SiiSili 
H(C) = H y (x, y) + iH x (x, y) as a function of ( = x + iy 
in the thin-strip limit, d/a — > 0, is expressed as 



n(0 = H a 



i 

2^ 



du 



K z {u) 
C-u' 



(1) 



where the magnetic field H a is applied parallel to the y 
axis, and K z (x) = j z (x,y)dy is the sheet current 

along the z axis. In Sees. IIVI and we show distribu- 
tions of magnetic field H y (x,Q) = Kc[H(x)] and current 
K z (x)/2 = tH x {x,±0) = =flm[H(x ± i0)]. The mul- 
tipole expansion of Eq. Q for ICI/ a ~~ * 00 is expressed 
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u 

7~2 



(2) 



where I z = dxK z (x) is the transport current along 



the z axis and m y = J_ a dx(—x)K z (x) is the magnetic 
moment in the y direction per unit length. In this paper 
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FIG. 1: Cross sections of strips with total width 2a and thick- 
ness d: (a) single strip without slit, (b) two strips (i.e., strip 
with a slit), and (c) three strips (i.e., strip with two slits). 



we consider the hysteretic relationship between m v and 
H a of strips carrying no net current (I z = 0) . 

The complex field for symmetrically arranged strips 
has the general form [e.g., Eqs. ©, pi|l. and 



C 2 



(3) 



where the strip edges are at x = ±a ra . Equation © may 
be expanded as 



no - if + 1| £ (o» 



Comparing Eq. (0J with Eq. J2J), we obtain a general 
expression for the magnetic moment per unit length [e.g., 
Eqs. 0, US), and 



(5) 




III. SINGLE STRIP WITHOUT SLITS 

A. Complex field for a single strip 

In this section we consider a single strip of width 2a, 
as shown in Fig. QJa). The complex field 7i(C) and the 
magnetic moment per unit length m y for a single strip 
are given by^ 

(6) 
(7) 

For convenience we introduce a± = a ± 5, where 8 
is a cutoff length on the order of the thickness d when 
the penetration depth A is small (A < d) or the two- 
dimensional screening length A = 2A 2 /c? when A is large 
(A > d); i.e., S ~ max(c?, A). 

B. Single strip in an ascending field 

Step (0) for < H a} < H 

When a magnetic field H a is increased after zero-field 
cooling, the superconducting strips are initially in the 
Meissner state. The magnetic field at the edges is less 
than the flux-entry field H s [i.e., H y (a+,0) < H s ], and 
the edge barriers near x ~ ±et prevent penetration of 
magnetic flux. In this step (0), the parameter in Eq. © 
is a — 0, and the magnetic moment is given by 



triy/'K = -H a a 2 . 



(8) 



Magnetic flux cannot penetrate into the strips so long as 
H y (a+,0) < H s for H a < Hq, but step (0) terminates 



(4) when H y (a + , 0) = H s at H a = H , 



Ho 

where the second equality is valid for 5/a <C 1. 

Step (i) for H <H a -t < H irT 



(9) 



When H a is small enough, such that the local mag- 
netic fields at the strip edges are all less than the flux- 
entry field 7? s! iMthe superconducting strips are in the 
Meissner state and no magnetic flux penetrates into the 
strips. References lol and ITol describe the linear reversible 
magnetic response of two strips in the Meissner state. On 
the other hand, when H a is sufficiently large to make the 
edge fields reach H s , magnetic flux penetrates into the 
strips, and the magnetic response becomes irreversible 
and hysteretic. In Sees. ITTTlVl we determine the parame- 
ters a n and the magnetic moment m y as functions of H a 
in ascending (H a ^) and descending (H a i) fields. 



When H a > Hq, magnetic flux nucleates at x ~ ±a 
and penetrates into the strip. A domelike distribution 
of magnetic flux exists for \x\ < ot\ and grows as H a 
increases. The parameter a — a.\ in Eq. JSJ) is determined 
by H v (a+,0) = H s , 

a\ = a 2 + -(H s /H a f{a\-a 2 ) 

~ a 2 ~2Sa[(H s /H a ) 2 -I]. (10) 

The magnetic moment is given by 

my/it = -H a (a 2 - a 2 ) ~ 2(H a - H 2 /H a )aS. (11) 
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Step (i) terminates when a% = a_ at H a = H lTl: , where 



strips, f_° K z dx = — K z dx ^ 0. General expres- 
sions for the complex field and the magnetic moment per 
unit length for two strips are given byi£ 



(12) 

4a V2 V ' 



Step (ii) for H a > Hi rl 

The domelike distribution of magnetic flux expands to 
include almost all of the strip, and the strip's magnetic 
response is reversible for H a > H 1TT . The edge field at 
H a = H irl is given by H y (a + , 0) ~ H s . 



no = H a 



'(C 2 -a 2 )(C 2 -/? 2 ) 



(C 2 -a 2 )(C 2 -6 2 ) ' 
= 7t H a (a 2 + b 2 - a 2 - P 2 ) 



(14) 
(15) 



In this section we use parameters a± — a±S and b± = 
b±S [where S ~ max(d, A)], and define the function 
hi^) — (x 2 — a 2 ) (a; 2 — b 2 ) for convenience. 



B. Two strips in an ascending field 



C. Single strip in a descending field 



Step (0) for < H a t < Ho 



Step (m) for < Hai < Hrx 

In a descending field, magnetic flux escapes from the 
strip, but a domelike distribution of magnetic flux re- 
mains in the strip. The detailed behavior of the field 
distributions and m y in descending fields depends upon 
the treatment of edges of a stripi^^ Here we adopt a 
simple model and put a — a_, which results in a dome- 
like distribution of magnetic flux for \x\ < a_ and a small 
but sharply peaked current density flowing in the vicin- 
ity of the edges, a_ < |a;| < a. The magnetic moment is 
given by 



my/-* = -H a (a 2 - a 2 _) 



-2H a aS, 



(13) 



which qualitatively agrees with that predicted by more 
detailed investigations!^^ 

Note that m y takes a very small negative value (i.e., 



< 



<C H s a ), because a single strip cannot trap 



any magnetic flux in a descending field. Geometric edge 
barriers in a single strip without slits do not prevent es- 
cape of magnetic flux. At H a = magnetic flux is en- 
tirely removed from a single strip, resulting in zero rema- 
nent magnetic moment, m y = at H a i = 0. For strips 
with slits, on the other hand, magnetic flux is trapped in 
the slits, and m y becomes positive in descending fields, 
as we show in Sees. II Vl andlvl 



IV. TWO STRIPS (STRIPS WITH A SLIT) 

A. Complex field for two strips 

In this section we consider two strips (i.e., a strip with 
a single slit) of total width 2a, as shown in Fig. ^b). 
The superconducting strips are at b < \x\ < a, and the 
slit is centered between the strips, \x\ < b. We assume 
that the two strips are connected at the ends (|z| — * oo) 
by superconducting shunts, such that an applied mag- 
netic field H a induces a circulating current in the two 



When a magnetic field H a is applied after zero-field 
cooling, the superconducting strips are initially in the 
Meissner state. The magnetic fields at the edges are less 
than the flux-entry field, < -H y (b^,0) < H y (a +1 0) < 
H s , and no magnetic flux penetrates into the strips. Fig- 
ure Et^a) shows distributions of the perpendicular field 
H y (x,0) and current K z {x)/2 = ^pH x (x, ±0) in step (0). 
The parameters in Eq. I|14|) are given by a = (3 = ctQ < 
b. Because the total magnetic flux in the slit is zero, 
H y (x,0)dx = 0, the parameter ao is determined as 



xl = a 2 [1 - E{b/a)/K(b/a)} , 



(16) 



where K(k) and E(k) are the complete elliptic integrals 
of the 1st and 2nd kind, respectively. The magnetic mo- 
ment is given by 



m y /TT = ~H a (a 2 + b 2 - 2«q). 



(17) 



Step (0) terminates when H y (a+, 0) = H s at H a = H , 
where 



Ho _ ^/Ua+) v/2oa(a 2 - 6 2 ) 



H, 



(18) 



Step (i) for H < H a f < #i 

During this step, magnetic flux nucleates at the outer 
edges (x = ±a), flows inward across the strips, and pen- 
etrates into the slit. Thus, we have H y (a+,0) = H s and 

H y (x, 0)dx > 0. The field distributions in step (i) are 
similar to those in Fig. |2Ia). The parameter a — (3 = a\ 
is determined by H y (a+,Q) = H s , 



a 2 + - (H t /H a )y/hfa 



- (H s /H a )^/2Sa(a2-b2). 



(19) 



The magnetic moment m y is given by 

my/n = -Hja 2 + b 2 - 2a\) 

~ H a (a 2 -b 2 ) ~2H sV / 25a(a 2 -b 2 ). (20) 
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FIG. 2: Field distributions of H y (x,0), shown as black lines, and K z (x)/2 — ^H x (x,±0), shown as gray lines, in two strips 
for which b/a = 0.5: (a) < H a] < Ho, (b) H a f = Hi, (c) Hi < H a i < H itI , (d) H 4 < H al < H iTr , (e) H al = Hi, and (f) 
H$ < H a i < Hi. The horizontal dot-dashed lines show the applied magnetic field H a . 



Step (i) terminates when a.\ = b+ at H a = Hi, where 

(21) 



Hi = VMa 
H s 



-6 



25a 



i 2 -b 2 ' 



The field distribution at H a = Hi is shown in Fig.|5fb). 



Step (ii) for Hi < H a i < H irl 

During this step, magnetic flux penetrates from x — 
±a, because H y (a + ,0) = H s . Domelike distributions of 
magnetic flux appear in the strips at 6+ < \x\ < a 2 
[Fig. 13(c)] , and grow as H a increases. Magnetic flux si- 
multaneously exits the superconducting strips at x = ±b 
and enters the slit. Current spikes occur near x ~ ±6 
in Fig. Efc) because a finite current K z ^ flows in the 
vicinity of the inner edges, b < \x\ < o+, whereas K z = 
for b+ < \x\ < 062. In our theory, such current spikes, 
which produce wiggles in the local magnetic field distri- 
bution, always occur where magnetic flux is exiting from 
a dome in a superconducting strip. 

The parameters are given by a = oli and /3 = &+, 



where a-i is determined by H y (a+, 0) = H s 

/2O+) 



<4 



H a ) a\-b\ 
~ a 2 -25a[(H s /H a ) 2 - 1]. 

The magnetic moment m y is given by 

m y /n - 



(22) 



2 , .2 2 .2 
O — a, — 0, 



H a {d T w 1 ! j 

2(H 2 /H a )aS + 2H a {a + b)S. (23) 



Step (ii) terminates when a2 = a_ at the irreversibility 
field -ff Q = i/i rr ~ H s /y/2, when domelike flux distribu- 
tions essentially fill the strips. The domes occupy the 
regions 0+ < |x| < a_. 



Step fm) /or _fT a > Hi Tr 

When H a = H m , domelike flux distributions fill the 
regions 6+ < \x\ < a_, and the magnetic response be- 
comes reversible. At H a — H m , the magnetic fields at 
the strip edges are H y (a+,0) ~ H y (b-,Q) ~ i/ s , and 
the parameters are given by a = o_ and [3 — b + . The 
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magnetic moment m y , which is due to currents flowing 
in the narrow regions near the edges (6 < |x| < b + and 
a_ < |x| < a), is given by 



ITly/lT 



-H m {a 2 + b 2 - 
-2H irr (a - b)S. 



i 2 - b\) 



(24) 



When H a > H 1TI , the magnitude of m y is reduced be- 
low that given in Eq. (|24[1 . because the current-carrying 
regions near the edges become narrower. However, a 
more detailed theory beyond the scope of the present ap- 
proach would be required to calculate m y for H a > Hi„. 



C. Two strips in a descending field 

Step (iv) for Hi < H al < H irl 

If the applied field H a has been above H m and now 
decreases through H m , magnetic flux is expelled from the 
slit and penetrates into the strips from the inner edge at 
x = ±b, because the inner-edge field H y (6_ , 0) is equal 
to the flux-entry field H s . Magnetic flux in the strips 
escapes from the outer edges at x = ±a. The domelike 
flux distributions at (3i < |x| < a_ shrink as H a decreases 
[Fig-Eld)]. Current spikes occur near x ~ ±a in Fig. Eld) 
because a finite current K z ^ flows in the vicinity of the 
outer edges, a_ < |x| < a, whereas K z = at 04, < |x| < 
a~. The parameters are given by a = a_ and /3 = 
where (3± is determined by H y (b-, 0) = H s , 



Pi 



= b 2 _ + 



HsY /a(6-) 
Hj a 2 _-b 2 _ 

b 2 + 28b[(H s /H a ) 2 - 1] 



(25) 



The magnetic moment m y is given by 

TFly 1 7T = 



-H a {a 2 + b 2 -a 2 _ -01) 
2(H 2 /H a )bS - 2H a (a + b)S. 



(26) 



Step (iv) terminates and the domelike field distribu- 
tions in the strips disappear when (3^ — a_ at H a = 
[Fig-Ete)], where 



h 4 _ ylM 

H s a 2 _ - b 2 _ 



2bS 



b 2 ' 



(27) 



The parameter a 
H y (b_,0) = H s , 



p 



as is determined by 



= b 2 



(H s /H a )^Mb- 



~ b 2 + (H s /H a )^25b(ai-b2). (28) 

Note that «5 is real (and > a_) for < H a < Hi, 
that as — * oo for H a — * +0, and that a§ is imaginary 
(i.e., a\ < 0) for H 5 < H a < 0. The magnetic moment 
m y is given by 



m y /ir 



H a (a 2 
-H a {a 2 



b 2 
b 2 ) 



2*1) 

f 2H sV / 2Sb(a 2 -b 2 ). (29) 



Even when H a < 0, positive magnetic flux is still trapped 
in the slit, and the magnetic moment is positive. 

Step (v) terminates when H y (a + ,0) = —H s at H a = 
H 5 < 0, 



-H 5 



1 



b 2 _ 



Vb 



25 



a 2 -b 2 ' 



(30) 



Step (vi) for -Hi < H al < H 5 

During this step, negative magnetic flux (i.e., flux lines 
aligned along the —y axis) penetrates into the strips from 
the outer edges (x = ±a), flows inward entirely across 
the strips, and annihilates at the inner edges (x = ±6), 
where H y (b-,0) > 0. The parameter a = f3 = a§ is de- 
termined by H y (a + , 0) = —H s , and is given by ae(H a ) = 
ai(— H a ), where a.\{H a ) is defined in Eq. I|19|l . Note that 
a6 is imaginary (i.e., a\ < 0) for H* < H a < H5, and 
that < ae < b + for —Hi < H a < H* , where Hi is 
defined in Eq. (|2*T)> and 



-H* 

~h7 



Vh(a] 



2(a 2 -b 2 )S 



(31) 



The magnetic field at the center of the slit, H(Q) — 
#3,(0,0) = H a a(3/ab, obeys H(0) > for H al > H*, 
and H(0) < for H al < H*. 



Behavior for H a i < —Hi 



Step (v) for H 5 < H al < Hi 

During this step, positive magnetic flux exits from the 
slit, penetrates into the strips from the inner edges at 
x = ±b, flows outward entirely across the strips, and an- 
nihilates at the outer edges (x = ±a), where H y (a + , 0) < 
0. The field distribution is shown in Fig.EtO- 



The magnetic response of two strips for H a [ < —Hi in 
descending fields is very similar to the response for iJ a | > 
Hi in ascending fields. The complex field H a ) 
and magnetic moment m y (H a ) in descending field can 
be determined from the ascending-field results with the 
help of the symmetries H((,H a i) = —H((,—H a ^) and 
m y (H al ) = —m y (-H a i), respectively. 
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D. Magnetization curves of two strips 

Figure |3 shows hysteretic m y versus H a curves for two 
strips. In ascending fields H a j, the sheet current K z (x) 
concentrates near the outer edges of the strips at x ~ ±a 
[Figs. Etb) and (c)], and therefore, even for large b/a, the 
m y for two strips is almost the same as that for a single 
strip except at low fields, H a /H s < 0.2. 

In descending magnetic fields H a ^ , we see in Fig. [3] a 
striking difference between m y > for two strips and 
m v w (i.e., < —m v <C H s a 2 ) for a single strip. The 
inner edges of the two strips at x = ±6 are responsible for 
the large positive magnetic moment. The edge barriers 
near x — ±6 prevent magnetic-flux penetration from the 
slit region |ac| < b into the superconducting region b < 
\x\ < a. Magnetic flux is therefore trapped in the slit 
as shown in Figs, ^d)-(f), and the magnetic moment 
becomes positive (ra y > 0). 

The remanent magnetic moment m rom is given by 
Eq. |29 with H a = 0, 



i/ 7T ~ 2H s y/25b(a 2 -b 2 ), 



(32) 



which is maximized when b/a — l/VS ~ 0.577. If the 
slit is made wider, m rcm increases as b/a increases for 
< b/a < 1/V3, because the slit can trap a larger 
amount of magnetic flux. If the slit is too wide, how- 
ever, m rem decreases for 1/V3 < b/a < 1, because the 
superconducting strips (i.e., the current-carrying region) 
become too narrow. A similar behavior occurs for the 
peak of the magnetic moment, m pca k, which occurs at 
H a = H 5 . Equations an d £23 yield 

m poak /7r ~ H s ( + Vbj v/2(a 2 - b 2 )6, (33) 

which is maximized when b/a~ 0.403. 



V. THREE STRIPS (STRIPS WITH TWO 
SLITS) 

A. Complex field for three strips 

In this section we consider three strips (i.e., strips with 
two slits) of total width 2a, as shown in Fig. ^c). The 
outer strips are at b < \x\ < a, the inner strip is at \x\ < c, 
and the slits are at c < \x\ < b. The strips are connected 
at the ends, \z\ — > oo, so that a circulating current flows 



in the outer strips, J_ a K z dx = — K z dx =/= 0. The 

inner strip carries no net current, f^ c K z dx = 0. The 
general expressions for the complex field and the mag- 
netic moment per unit length for three strips are 



no = H ai 



'(C 2 -a 2 )(C 2 -/3 2 )(C 2 -7 2 ) 



(C 2 -a 2 )(C 2 



(34) 



-nH n 



b 2 + c 2 



b 2 )(( 2 -c 2 ) ' 
- a 2 - (3 2 - 7 2 ) . (35) 
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FIG. 3: Hysteretic behavior of the magnetic moment m y of 
two strips as a function of the applied magnetic field H a . 
Magnetization curves of a single strip (dashed lines) and two 
strips (solid lines) for b/a — 0.1, 0.577, and 0.9. Thickness of 
the strips is d/a — 0.01. 



In this section we use parameters, a± — a ± 5, b± = 
b ± S, and c± — c ± S, where 5 ~ max(<i, A). We also 
define a function f^{x) = (x 2 — a 2 )(x 2 — b 2 ){x 2 — c 2 ) for 
convenience. 



B. Three strips in an ascending field 

Step (0) for < H a] < Ho 

When a magnetic field H a is applied after zero-field 
cooling [step (0)], the superconducting strips are initially 
in the Meissner state. Figure^a) shows the distributions 
of H y (x,0) and K z (x)/2 = ^fH x (x,±0). The parameters 
in Eq. I|34ll are a = /3 = ao and 7 = 0, where c < ctQ < 
b. Because the net magnetic flux in the slits is zero, 



J H y {x, 0)dx — 0, we find that 



a 2 = a 2 - (a 2 - c 2 ) E{k) 



k = 



b 2 



K(k)' V a- 

The magnetic moment is 

m y /iT = -H a (a 2 + b 2 + c 2 - 2al). 



(36) 



(37) 



Step (0) terminates when H y (a + , 0) = H s at H a = Ho, 
where 



Hp _ \/h{a+) 
H, 



a+(a 2 



{)) 



^J25a{a 2 - b 2 ) {a 2 ~ c 2 ) 
a{a 2 — q.q) 



(38) 



Step (i) for H < H a f < (H 1A or Hm) 

In this step, magnetic flux nucleates at the outer 
edges, x = ±a, flows inward across the outer strips, and 
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FIG. 4: Field distributions of H y {x,Q), shown as black lines, and K z (x)/2 = ^H x (x, ±0), shown as gray lines, in three strips 
for which b/a = 0.7 and c/o = 0.4: (a) < H a} < H , (b) H 1A < H a f < H 2A , (c) # 2 a < #at < #irr, (d) H s < H ai < H ilT , 
(e) H§ < H a i < an d (f) ^8 < H a i < (magnified distribution in the inset). The horizontal dot-dashed lines show the 
applied magnetic field H a - 



penetrates into the slits; that is, H y (a+,Q) = H s and 
J c H y {x, 0)dx > 0. However, no magnetic flux penetrates 
into the inner strip because H y (c+, 0) < H s . The pa- 
rameters are given by a = (3 = a\ and 7 = 0, where 
c < ct\ < b + . The value of ct\ is determined from 
H y (a+,0) = H S , 



- (H s /H a )^M^)/a + 



~ a 2 - (H s /H a )^2(a2 - fo 2 )(a 2 - c^)5/a. (39) 

The magnetic moment m y is given by 

my/n = -H a (a 2 + b 2 +c 2 -2a\) 
~ +H a (a 2 -b 2 - c 2 ) 



-2H s yj2{a? - 6 2 )(a 2 - c 2 )6/a. (40) 

The sign of dm y /dH a = n(a 2 — b 2 — c 2 ) can be either 
positive or negative, depending upon the relative widths 
of the strips. 

Step (i) terminates when either a\ — b + at H a = Hia 
or H y (c+ ,0) = H s at H a = Hib ■ The characteristic 
fields -Hia and Hib are given by 



H 



1A 



~h7 



a + (a 2 + 
1 



l2{a 2 -c 2 )S 
a(a 2 - b 2 ) ' 



(41) 



a+ 



c+ 



8 



28 



a 2 -b 2 



b 2 -c 2 



(42) 



When b 2 < a 2 — ac + c 2 , step (i) terminates at H a = 
Hia < Hib , whereas when b 2 > a 2 — ac + c 2 , step (i) 
terminates at H a = Hib < -Hia- 



Step (ii-A) for H\a < H a f < H 2 a and b 2 < a 2 — ac + c 2 

During this step, domelike distributions of magnetic 
flux are present at 6+ < |x| < ot2A in the outer strips, 
whereas no magnetic flux penetrates into the inner strip 
[Fig. Eft))]. The parameters are given by a — «2A, 
(3 = b + , and 7 = 0, where b + < a.2A < and a2A 
is determined by H y {a+, 0) = H s , 



"2A 



Ha+) 



_(a%-b%) 



~ a 2 -2(H S /H a ) 2 (a 2 -c 2 )8/a. 
The magnetic moment m y is given by 



(43) 



m y /7T = -H a (a 2 +b 2 + c 2 



l 2A ~ u + 



~ -H a c 2 -2(H 2 /H a )(a 2 - c 2 )S/a. (44) 

Step (ii-A) terminates when H y (c + ,0) = H s at H a = 
H 2 a, where 



Ho 



1 



fs(c+) 



ai(a 2 + -b 2 + ) c 2 + (b 2 + -c 2 + ) 



y/2(a + c)8/ac. 



(45) 



Step (ii-B) terminates when a2B = b + at H a = H2B, 
where 



H2B 



1 



/3(a+) /3(c+) 



(al-bl) 2 {b\-clf 



'2(a + c)(6 2 -ac)6 
(a 2 -b 2 ){b 2 -c 2 )' 



(48) 



Step (m) for (H 2 a or H 2 b) < H a ^ < H irl 

During this step, domelike distributions of magnetic 
flux arc present both for |x| < 73 in the inner strip and 
b+ < \x\ < 0:3 in the outer strips [Fig.^Jc)]. The param- 
eters are given by a = a.3, — b + , and 7 = 73, where 
< 73 < c_, b + < a 3 < a-, and 03 and 73 are de- 
termined by the coupled equations H y (a + ,0) = H s and 
H y (c + ,0) = H s , which yield 

H s /H a = ^{a 2 + -a 2 ){a 2 + - 1 l)(a 2 + -b 2 + )/h{a + ) 



xj(a 2 - c\){c 2 + - j 2 )(b 2 + - c\)/h{c + ). 



(49) 



The magnetic moment is given by 
m y 1 7r = -H a (a 2 + b 2 + c 2 - aj 



~ 2H a (a + b + c:)5- 2(H 2 /H a )(a + c)8. (50) 

Step (iii) terminates when 03 ~ a_ and 73 ~ c_ at 
H a = H m ~ H s /^/2. 



Step (iv) for H a > H ir 



Step (ii-B) for Hib < H a -\ < H2B and b 2 > a 2 — ac + c 2 

In this step, a domelike distribution of magnetic flux 
is present for \x\ < 72B in the inner strip, whereas no 
magnetic flux is present in the outer strips. The param- 
eters are given by a = (3 = ot2B and 7 = 72B, where 
< 72B < c_, c < a2B < b, and a2B and 72B are de- 
termined by the coupled equations H y (a+ ,0) = H s and 
H y (c + ,0) — H s , which yield 



H s /H a = {a 2 + - a 2 B)^(a 2 + - 1 2 B)lh{a + ) 

= («2B-4)\/( c +- 72 2 b)//3(c + ). (46) 
The magnetic moment m y is given by 

m y /ir = -H a (a 2 + b 2 + c 2 - 2a 2 B - 7 2 2 B ), (47) 

where a2B and 72B must be determined numerically from 
Eq. P>|). 



The magnetic response is reversible for H a > H m . At 
H a = H lrl the parameters are a = a_, (3 = b + , and 7 = 
c_, and the magnetic fields at the edges are H y (a+, 0) ~ 
H y (b-,0) ~ H y (c+,0) ~ H s . The magnetic moment at 
H a = H m is given by 



nXyj-K = — H- 1XT {a 2 + b 2 + c 2 
~ -2H iTt {a - b + c)5. 



c 2 -) 



(51) 



C. Three strips in a descending field 

Step (v) for H 5 < H al < H irl 

During this step, domelike distributions of magnetic 
flux are present for (3$ < \x\ < a_ in the outer strips and 
for \x\ < c_ in the inner strip [Fig.Hfd)]. These domelike 
flux distributions shrink as H a decreases. Magnetic flux 
escapes from the inner strip and penetrates into the slits. 
In turn, magnetic flux exits from the slits, penetrates into 
the outer strips, flows outward along the outer strips, 
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and finally escapes from the strips at the outer edges, 
x = ±a. The parameters are given by a — a_, (3 — (3$, 
and 7 = c_, where is determined by H y (b-, 0) = i? s . 



# = + ( |± 



/s(6-) 



Hj (a 2 _ - 6 2 _)(6 2 _ - <?_) 
~ 6 2 + 2^[(ff s / J ff a ) 2 -l]. (52) 

The magnetic moment m„ is given by 

m^/TT = -i? a (a 2 + & 2 + c 2 -a 2 _-/3|-c 2 _) 

~ ~2H a (a + b + c)6 + 2{H 2 /H a )b6. (53) 

Step (v) terminates when /3g = a_ at _ff a = where 

#5 



_ _ 1 My 

# s a 2 - 6 2 _ V fe 2 _ - c 2 . 



265 



& 2, 



(54) 



Siep fw,) /or H 6 < H al < H s 

Throughout this step, no magnetic flux remains in the 
outer strips, but a domelike distribution of magnetic flux 
is still present for \x\ < c_ in the inner strip [Fig. life)]. 
Magnetic flux at |x| < b penetrates into the outer strips 
from x = ±6, flows outward entirely across the outer 
strips, and annihilates at x = ±a, because H y (a+, 0) < 0. 
The parameters are given by a = (3 = a§ and 7 = c_, 
where ae is determined by H y (b^, 0) = H s , 



fs(b-) 



H a yb 2 _- C *_ 



b 2 + (H s /H a )^/25b(a 2 -b 2 ). 



(55) 



Note that is real (and > a) for < H a < H 5 , that 
a — > +00 for H a — > +0, and that a 6 is imaginary (i.e., 
ag < 0) for Hq < H a < 0. The magnetic moment m y is 
given by 

m y /7r = -H a (a 2 + b 2 + c 2 - 2a% - c 2 _) 

~ -i/ a (a 2 - 6 2 ) + 2H sy /28b(a 2 -b 2 ). (56) 

In the thin-strip limit d/a — > 0, the remanent magnetic 
moment m r em 5 given by Eq. 1)56(1 with H a — 0, is the 
same as for two strips, Eq. I|32|) . 

Step (vi) terminates when H y (a+,0) = —H s at H a = 
H 6 < 0, where 



-g 6 



/3(a+) 



fs(b- 



2(5 



6 2 ' 



(57) 



Siep (wij /or #7 < H ai < He < 

During this step (not shown in Fig. 4), the magnetic 
field at the outermost edges is equal to the negative flux- 
entry field, H y (a + ,0) = —H g . Negative magnetic flux 
penetrates into the outer strips at x = ±a and flows 
entirely across the outer strips into the slits, resulting in a 
reduction of magnetic flux in both the slits and the center 
strip. The parameter 7 = c_ is fixed, and a = (3 = a 7 
is imaginary (i.e., a 2 < 0), where aj is determined by 
H y (a+,0) = -H s , 



a 7 — a. 



H a \l a 2 -c 2 



~ a 2 + (H s /H a )^25a(a 2 -b 2 ). (58) 
The magnetic moment is given by 

m y /Tr = -H a (a 2 + b 2 + c 2 - 2a 2 7 - c 2 _) 



~ +H a (a 2 - b 2 ) + 2H s ^/2Sa(a 2 - b 2 ). (59) 

Step (vii) terminates when aj = at H a = Hj < 0, 
where 



-H r 



1 / 3 (o+) 



2(a 2 - b 2 )8 



(60) 



The magnetic field at the center, H(0) — H y (0,0), is 
positive for H a > Hj, but becomes zero at H a = H7. 



Step (viii) for H 8 < H al < H 7 < 

During this step, domelike distributions of magnetic 
flux are present at /3g < \x\ < c_ in the inner strip, 
whereas no magnetic flux is present in the inner re- 
gion, |x| < (3s [Fig. Elf)]- The parameters are given by 
a = 0, f3 = 13%, and 7 = c_, where f}% is determined by 
H y {a+,Q) = -H s , 



4 



H a J a 2 + (a 2 + -c 2 _) 



~ a 2 -2{H s /H a ) 2 {a 2 -b 2 )S/a. (61) 

The magnetic moment is given by 

m y /ir = -H a (a 2 + b 2 + c 2 - (3 2 - c 2 _) 

~ -H a b 2 ~2(H 2 /H a )(a 2 -b 2 )S/a. (62) 

Step (viii) terminates when /3g — c_ at H a = Hs < 0, 
where 



-H 8 
H, 



VJW^) 2(a 2 -b 2 )S 



a + (a + — c_) y a(a 2 — c 2 ) 



(63) 



At H a = H s , no magnetic flux is present in either the 
inner or outer strips. 
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Step (ix) for -Hi < H al < Hs < 

During this step, no magnetic flux is present in the 
strips, and the parameters are given by a = and (3 = 
7 = /3q. The parameter /J g , determined by H y (a + ,0) = 
—H s , is given by /3g(H a ) — o>x(—H a ), where a\ is defined 
in Eq. 

Behavior for H a i < —Hi < 

The magnetic response of two strips for H a i < —Hi in 
a descending field is very similar to that for H a ] > Hi 
in an ascending field. The complex field and magnetic 
moment as functions of the applied field can be deter- 
mined with the help of the symmetries H(C, H a i) — 
—H((,—Hat) and m y (H a i) = —m y (—H a t), respectively. 

D. Magnetization curves of three strips 

Figure shows hysteretic m y -H a curves for three 
strips. In ascending magnetic fields at H a i = 
(i?2A or i?2B), the m y -H a curves of three strips have ad- 
ditional kinks [see arrow in Fig.^a)] or peaks [see arrow 
in Fig. EJb)] . Such kink or peak structures do not ap- 
pear in the m y -H a curves for a single strip or two strips. 
They arise in three strips because the the edge barrier at 
the edges of the center strip (x ~ ±c) impede the entry 
of magnetic flux. In descending fields H a ^, on the other 
hand, the m y -H a curves for three strips are almost the 
same as those for two strips; the similarity occurs because 
the edges of the center strip do not impede the exit of 
magnetic flux. 

In three strips, the edge barriers near x ~ ±a and ±c 
are effective in preventing the entry of magnetic flux in 
H a i > but are not effective in stopping the exit of magnetic 
flux in H a i ■ On the other hand, while the edges at x ~ ±b 
do not impede entering magnetic flux in H a -\, the edge 
barriers there are responsible for impeding flux exit in 

Hal- 

VI. CONCLUSION 

We investigated field distributions and the magnetic 
moment m y of bulk-pinning-frcc strips with slits in ap- 
plied magnetic fields H a , and we studied these in detail 
for increasing fields H a -\ and decreasing fields H a i- 

For two strips in decreasing fields H a i, the edge bar- 
riers near the inner edges at x ~ ±6 impede the exit of 
magnetic flux from the slit at |a;| < b into the strip at 
b < \x\ < a. The trapping of magnetic flux in the slit 
in H a ^ results in a positive remanent magnetic moment, 
m y > at Hai = 0> rather than a zero remanent mo- 
ment, which occurs for a single strip without slits. For 
two strips, the remanent magnetic moment at H a i = 
is maximized when b/a — l/v3. 




-0.2 0.2 0.4 0.6 0.8 

FIG. 5: Hysteretic behavior of the magnetic moment m y of 
three strips as a function of the applied magnetic field H a . 
(a) Magnetization curves of a single strip (dashed lines), two 
strips (gray solid lines) with b/a = 0.7, and three strips (black 
solid lines) with b/a = 0.7 and c/a — 0.5. (b) Magnetization 
curves as (a) but for b/a — 0.9 and c/a = 0.8. Thickness of 
the strips is d/a = 0.01. 



For three strips (i.e., strips at \x\ < c and b < \x\ < 
a) the edge barriers at x ~ ±b in H a i are effective in 
trapping magnetic flux as in two strips. In the thin-strip 
limit d/a — > 0, the remanent magnetic moment for three 
strips at H a i = is the same as for two strips. In H a i, 
the edge barriers at x ~ ±c impede the penetration of 
magnetic flux into the inner strip. As a consequence, the 
m y of three strips in H a ] exhibits kink or peak structures 
that are not present for two strips. 

The above arguments for three strips can be extended 
to an arbitrary number of strips. Consider a symmetric 
array of TV coplanar strips of total width 2ao (N — 1 
slits), in which the strip edges are at x = ±o„, where 
cij\[ —i < cln —2 < ■ ■ ■ < ai < do- Superconducting strips 
occupy the regions a n+ i < |x| < a n with n even, and 
slits are at a n+ i < \x\ < a n with n odd. The region 
spanning the centerline, |a;| < ajy-i corresponds to a 
superconducting strip for odd TV and a slit for even N. 
The even edges at x = ±a n for n = 0, 2, 4, • • ■ impede 
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flux entry in increasing applied fields H a ^, and the odd 
edges at x = ±a„ for n — 1, 3, 5, • • • impede flux exit in 
decreasing fields H a i- 

By considering narrow slit widths, one may further ex- 
tend the above approach to model grain-boundary pin- 
ning in thin films. The resulting model would be the 
thin-film analog of the bulk-pinning model discussed in 
Ref.El 
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